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Abstract 

We have performed simulations to study how increasing humidity affects the structure of Nafion- 
like ionomers under conditions of low sulfonate concentration and low humidity. At the onset of 
membrane hydration, the clusters split into smaller parts. These subsequently swell, but then 
maintain constant the number of sulfonates per cluster. We find that the distribution of water 
in low-sulfonate membranes depends strongly on the sulfonate concentration. For a relatively low 
sulfonate concentration, nearly all the side-chain terminal groups are within cluster formations, 
and the average water loading per cluster matches the water content of membrane. However, for 
a relatively higher sulfonate concentration the water-to-sulfonate ratio becomes non-uniform. The 
clusters become wetter, while the inter-cluster bridges become drier. We note the formation of 
unusual shells of water-rich material that surround the sulfonate clusters. 
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PACS: 61.41. +e; 82.47.Nj; 64.75. St; 82.47.Gh 



I. INTRODUCTION 

An interest in ionomers, i.e. ion-containing polymers, began fifty years ago with the 
development of organic ion-exchange resins jl|. The properties of these materials are com- 
pletely different from those of other polymers as a consequence of the ionization of the ionic 
groups in polar solvents. The specificity of the interaction between the ion, the solvent, and 
the polymer makes it possible for these materials to be used as perm-selective membranes, 
thermoplastics or films for micro-encapsulation and coating A significant interest in 
ionomer materials also stems from their growing application as a polymer electrolyte mem- 
brane (PEM) or a proton exchange membrane in fuel cell technology 

DuPont de Nemours was the first manufacturer in the early 1960s to develop a perfiuoro- 

nn. 

sulfonic membrane commercially [7|, 18|. This membrane, which was named Nafion®, consists 
of a polytetrafluoroethylene hydrophobic backbone to which perfiuorovinyl ether pendant 
side chains are attached at more or less equally spaced intervals. The pendant chains are 
terminated by sulfonic head groups SO3H, and these are responsible for the large variety 
of microstructures in which the ionomer can be assembled. When exposed to humidity, the 
membrane takes up large amounts of water, leading to the dissociation of the acid groups 
SOsH^SOs + H+ and to the formation of a nanophase-separated network of aqueous (hy- 
drophilic) clusters and hydrophobic polymer. According to the cluster morphology model 
of Hsu and Gierke j^, spherical clusters are uniformly distributed throughout the material 



and are interconnected by channels 10, UJl- Subsequent cluster-based models, such as the 
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Mauritz-Hopfinger model [12|, 113||, the Yeager three-phase model the Eisenberg model 
of clusters of hydrocarbon ionomers jl5|, and the Litt model of a lamellar morphology for 
sulfonate domains [l^, have tried to quantify the cluster radius and spacing as a function of 
the polymer equivalent weight and the hydration level. Other structural models were pro- 



posed to describe the mem 

n 

as a continuous process [17 



Drane swelling process from a dry state to a colloidal suspension 



There is still ongoing debate about which one of the proposed models is more suitable 
and effective in representing the ionomer's conductivity through its nanophase separated 
network of hydrophilic regions. The issue is complicated by the fact that experimental 
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studies show t 
preparation [8|, 



le ionomer structure to depend on the pretreatment methods used in its 
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23|. The membrane pretreatment serves to reduce the remnant 



anisotropy in the morphology of extruded membranes, and to clean a solvent-cast membrane 



from impurities 
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2J|. Rigorously speaking, the question of how the pretreatment steps. 



such as swelling and/or boiling in solvents, annealing, rinsing in water, drying in vacuum/air, 
and the order of these steps, affect the membrane morphology, is not yet answered. Most 
of the pretreatment protocols have the ultimate goal of improving the water uptake of the 
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membrane [20, l25|. For example, in Ref. [20|| it has been shown that the water uptake of a 
dry membrane depends on how it was dried from its swollen state at elevated temperatures. 
If it was first cooled and then dried, then the membrane keeps its swollen volume. But if it 
was first dried and then cooled, the membrane shrinks in volume during the drying process. 
As a result, the outcome of the first protocol is a membrane that takes up a desirably large 
amount of water, and thus has a better proton conductivity. 

The water solvation of a PEM, which is necessary for its efficient operation, reduces its 
working temperature range: the membrane will not be exploitable at freezing and boiling 
water temperatures. A possible way to overcome this limitation is the development of new 
membranes that can operate at the low wetting conditions where A, which is the number 
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hydration there 
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28|. In low- 



of water molecules per sulfonate group, is less than five. In the case o: 
are 5 water molecules in the primary hydration shell of a sulfonate 
humidity membranes, the protons diffuse along narrow pathways near the SO3 terminals of 
side chains, and two confiicting effects come into play. On the one hand, the proximity of 
negatively charged sulfonates considerably suppresses the mobility of the protons. On the 
other hand, when the separation distance between sulfonates is small, the activation energy 
for proton hopping between adjacent end- groups becomes comparable with the activation 



energy in the bulk water 



29|, making the net result unclear. The proton mobility in low- 



humidity membranes can be also elevated by adding fiexibility to the sidechains, and by 
modifying the network structure of sidechain clusters. 

Despite the fact that various models have emerged to explain the properties of hydrated 
Nafion membranes, a systematic study of how the molar concentration r] of sulfonate head 
groups and the solvent content parameter A affect the network structure of sulfonates, and 
particularly the swelling of single clusters has not yet appeared. This absence is important 
for understanding proton transport and the onset of percolation in low-humidity membranes. 
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and forms the motivation for this study. 

Here we perform simulations to investigate the dependence of the cluster swelling on the 
hydration level A and the sulfonate molar concentration 7] of the membrane by employing dif- 
ferent sidechain architecture models. We restrict ourselves to the case of ionomers for which 
there is no bulk water inside the sulfonate clusters, and in which the sulfonate concentration 
is considerably below the percolation limit for the head groups. It is expected that in these 
low-humidity and low-sulfonate membranes no overlapping between sulfonate clusters takes 
place. In order to distinguish a sulfonate cluster from a water cluster, which is necessary 
in the interpretation of our simulation results, for the former we adopt the term "sulfonate 



multiplet", first introduced by Eisenberg in Ref. [15|| to describe the primary aggregates of 
sulfonates. We show that at the onset of swelling, which is defined as the transition from a 
dry multiplet into a wet multiplet with dissociated protons, the multiplets split into smaller 
parts. The solvation of these resultant multiplets is analyzed for different hydration levels 
and sulfonate concentrations. In particular, we will demonstrate the formation of water 
shells around the sulfonate multiplets. 

The paper is organized as follows. In Section II we briefiy discuss the benefits of using 
coarse-grained models (as opposed to all-atomistic approaches), and describe the coarse- 
grained model and system parameters employed here. The simulation details are outlined 
in Section III. Results on multiplet formation in dry and solvated membranes, water shells 
around multiplets, ionomer deformation, and proton diffusion are discussed in Section IV. 
We conclude in Section V. 



II. COARSE-GRAINED SYSTEM PARAMETERS 



Despite the rich variety of experimental findings and theoretical predictions for the or- 
dered morphology in PEM materials, numerical experiments have so far had little success 
in finding any clear picture of cluster formation in hydrated membranes. The main reason 
for this is the fact that individual ionic clusters are about 2-5 nm in size, and this is usually 
comparable to, or even larger than, the system sizes affordable in all-atomistic modeling. 
As a result, atomistic simulations, which are quite helpful for understanding the simple pore 
physics and small ionomer molecular conformations, are not able to capture the distribu- 
tion of sulfonate clusters in the hydrophobic matrix. However, as already outlined in the 
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introduction, a knowledge of this distribution is crucial for the determination of the ionomer 
connectivity and the proton conductivity of the PEM material. 

Fortunately polymers show a large degree of universality in their static and dynamic 
behavior. The universal scaling properties of the ionomer as a function of chain length, sul- 
fonate density, and membrane composition can be most efficiently studied via coarse-grained 
molecular models. One of the most commonly employed systems is a bead-spring model, 
where each bead represents a segment of a realistic chain. Wescott et al [s^ and Vishnyakov 
et al [3l| have simulated large ionomer systems using coarse-grained approaches in which 
an entire sidechain was represented by a nanometer-size hydrophilic blob. Their simula- 
tions report irregularly shaped hydrophilic clusters embedded into the polymeric matrix of 
backbone chains. While such gross coarse-graining is computationally convenient, it is not 
possible to draw firm conclusions regarding proton diffusion from the conformational results 
obtained for the polymer. It is therefore necessary to limit the coarse-graining approach to 
the level at which the sulfonic acid groups of the polymer can be explicitly treated, as these 
groups contain the essential membrane-specific interaction sites relevant to absorbed water 
and conducting protons. 

In our 'united atom' approximation for Nafion, the ether oxygens and sulfur atoms are 
treated individually, while the fluorocarbon groups are consolidated as a single particle, as 



are the three oxygens of the sulfonate 
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3j|. The fluorocarbon groups, the sulfonate 



oxygens, and the sulfur atoms are modeled as single Lennard- Jones (LJ) particles with a 
diameter a= 0.35 nm. The protons carry the full formal charge of = +e, the sulfur 
atoms have Qs = +l.le, and the combined triplet of oxygen atoms carries Q03 = — 2.1e. 
The partial charges of the ether oxygens and the fluorocarbon LJ particles are set to zero. 
Depending on whether the membrane is dry or hydrated, two different representations have 
been used for the sulfonate head groups. For dry membranes, we implement an attached- 
proton model, also called a dipole model for head groups, which was extensively analyzed 



m our previous paper 



32]. Though the attached-proton model does not allow for proton 



diffusion, it is considered as a good starting point for a step-by-step exploration of nanophase 
morphology in PEM materials. For the hydrated ionomer we assume a detached-proton 



model I33 



35|, where the protons diffuse freely in the system, where they interact with 



ionized head groups and water molecules. 

The configurational part of the coarse-grained Hamiltonian for the attached and detached 
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proton models is a combination of Coulomb interactions, non-bonded, and bonded interac- 
tions between all the ionomer components: 

fAotal = UlJ + Uq + [/bond + [Wangle + [^dihedral • (1) 

Here f/totai, Ulj, Uq, f/bond, f^angie and f/dihcdrai are the total potential energy and its Lennard- 
Jones, electrostatic, bond-stretching (bond-length term), angle bending (bond-angle term) 
and dihedral angle components, respectively: 

UlAt) = AsLj E ({^/r^,Y' - a{a/n,f) , (2) 

Uq = T.^^ (3) 

U^.oUr) = \ E h{r~r,)\ (4) 
all bonds 

u.nAd) = \ E ke{e-e,)\ (5) 

all angles 

[/dihedral(a) = ^ E K {I - d cof^{?,a)) . (6) 
all dihedrals 

In Eq.(l2]) the LJ interaction coefficient Slj in units of ksT was chosen to be 0.33. The 

parameter a in the LJ term is 1 for hydrophobic-hydrophobic (HH) interactions, and 0.5 for 

hydrophobic-hydrophilic (HP) interactions and hydrophilic-hydrophilic (PP) interactions. 

In the latter case only the repulsive part of the LJ potential (a shifted LJ potential for 

r < 1.1224cr) has been considered. In Eq.(l3]), Qi and Qj are the electrostatic charges 

of the two interacting elements, which can be sulfur atoms, oxygen triplets, protons, or 

the hydrogens or oxygen of the water molecules, and e is the dielectric constant of the 

ionomer. In Eqs.Q-Q the following force-field parameters have been used: the equilibrium 

bending angle 6q = 110", the equilibrium bond length tq = 0.44o", the bending force constant 

ke = 120—^4—1 and the stretching force constant kb = 7 x 10^ — tt^^?. The dihedral 
mol deg mol(nm) 

angle parameters were = — 1 (+1) and ka = lO.S/c^T {ka = S.TksT) for the backbone 
(sidechain) segments. 

The dielectric properties of the coarse-grained material are represented by a distance- 
dependent dielectric function, 

e{r) = l + eB{l-r/a)yy{l + ir/ar), (7) 



where is the bulk dielectric constants of the ionomer. Usually a uniform permittivity 
e = 1, or equivalently = in Eq.(l7l), is accepted in ab initio quantum-mechanical simula- 
tions, where all the ionomer atoms are explicitly taken into account. Since the coarse-grained 
approach neglects the atomistic structure of the ionomer monomers, additional approxima- 
tions for the dielectric permittivity have to be made to account for the polarization effects of 
the ionomer monomers as a response to the strong electrostatic fields of the sulfonate groups 
and protons. To be accurate, e(r) should depend upon the atom types and the absolute val- 
ues of all the explicit coordinates. The problem, of course, is that the specific form of e(r) 
is not known. For the bulk dielectric constant we use €3=8, which is appropriate to the 



19|| and differential 



dielectric permittivity of Nation as measured in high-frequency studies 
scanning calorimetry js^, and from first-principle calculations 37 1. 

Taking into account the fact that the sulfonic acid tips of sidechains are hydrophilic, and 



the remaining part of sidechains, as well as the backbone polymer, are hydrophobic [38|, 139|, 
we use the following notation to describe the polymer architecture: niH + n2P for sidechains 
and n^H for the backbone segments. Here ni is the number of hydrophobic monomers per 
sidechain, n2 is the number of hydrophilic monomers per sidechain, and is the number 
of backbone monomers between two adjacent sidechains. The total number of sidechain 
monomers per pendant chain is ni + ^2 in the detached proton model, and ni + ^2 + 1 in 
the attached-proton model. The key variables that describe our model ionomer system are 
listed in Table [B 



III. SIMULATION DETAILS 



Extensive coarse-grained molecular dynamics simulations were performed to investigate 
the swelling properties of sulfonate multiplets at four different solvation parameters A and 
two distinct sulfonate molar concentrations rj. The parameter rj is defined asr] = (Ns/Nq) V, 
where Ns is the number of sidechains in the volume V = of the simulation cell, and Nq 
is Avogadro's number. Whereas in experimental studies the parameters r] and A are coupled 
to each other [4^, in numerical simulations both quantities can be changed independently. 
A series of simulation runs are summarized in Table HTl The molar concentrations ?7i=0.8 
mol/1 and 772=!. 5 mol/1 correspond to the backbone segment lengths ^3=50 and 11^=20 
respectively. Varying the parameter r], i.e. varying the chain volume per SO3 group, is 
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in some ways equivalent to simulating materials with different equivalent weights {4^. For 
convenience, we will refer to the membrane with sulfonic molar concentration rj as 'membrane 
77'. In most simulations of Nafion-like ionomers the parameter is usually varied between 
14 and 18. Shorter segments with ^3=10 have been considered in the atomistic simulations 



of Ref. 



4l| . and longer segments ^3=30 in the coarse-grained approaches of Ref. [31|. In the 



latter case the nearest-neighbor distance between the sulfonate multiplets is large. Thus our 
choice of a larger makes possible the investigation of the solvation properties of single 
multiplets in slightly hydrated membranes. 

There were = (^1 + ^2 + ^^3) x Ns polymer monomers in the simulation box of length L. 
All simulations were carried out for both the rii = 7 and ^2 = 2 sidechain architectures. The 
number of sidechains Ns was 500 for membrane rji and 1000 for membrane ri2. The negative 
charges of A^^ sulfonate groups were compensated by Ns positive protons to guarantee an 
overall charge neutrality in the system. Simulations with explicit water include an additional 

nn 

3 X A X Ns water charges, as we are using the SPC solvent model [42, 143||. The box size 
L was systematically increased from L = 30a to L = 32.5a when the water content A was 
increased from A=0 to A=5 in order to keep the density of the hydrated membrane constant. 

One of the main challenges in generic ionomer simulations is the fact that the ionomer 
molecule is quite stiff at ambient temperatures and low humidity conditions. In experimen- 
tal studies a fast ionomer equilibration is usually achieved through different pretreatments 
protocols, such as a soaking in a solvent or high-temperature annealing. These steps im- 
prove the sidechain kinetics and decrease the barrier between the trapped metastable states 
and the low-lying states at the global minimum in free energy. Overall, a full equilibration, 
even after these pretreatment steps, takes hours or days, a time span that is far beyond 
the feasible simulation times of several nanoseconds in typical molecular dynamics runs. To 

n 

overcome this obstacle we implemented the following artificial steps |44l |: 



a) the sidechains were temporarily detached 



has been successfully applied in Refs. 



35 



rom the backbone skeleton, a technique that 



45|, 



46, 



b) the skeleton was cut into smaller segments of length n^. 

The resulting 'fragmented' ionomer reaches the equilibrium state very fast because of the 
increased diffusive movement of its segments. Typical MD runs of 500 ps duration in the 
NVT ensemble were enough to fully equilibrate the simulated system. The system temper- 
ature T was controlled by coupling the ionomer to a Langevin thermostat with a friction 
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coefficient 7 = 10ps~^ and a Gaussian white-noise force of strength GkBTj. The equations 
of motion were integrated using the velocity Verlet algorithm with a time step of 0.25 fs. We 
also imposed standard periodic boundary conditions to our system, thus filling space with 
translational replications of a fundamental cell. Long-range electrostatic interactions were 
treated using the Lekner summation algorithm 4^. 

In the next stage of the simulations, the ionomer segments were reassembled back into 
a branched chain characterizing the original Nafion-like ionomer. This was achieved by a 
simultaneous introduction of bonds and angular constraints between the ends of backbone 
segments unifying them into a single and long backbone chain. Similar bond and angular 
constraints were introduced between the fiuorocarbon tail monomers of detached sidechains 
and the median section monomers of backbone segments. To avoid the formation of unphys- 
ical star-like branches only a single occupancy of backbone attachment sites was permitted. 
The simulations were then resumed for another few hundred picoseconds until a new equilib- 
rium state was reached. Then the statistically averaged quantities of interest were gathered 
during the next 3ns-5ns of the long production runs. 

IV. SIMULATION RESULTS 

A. Multiplet formations in dry and solvated membranes 

A typical snapshot of a hydrated membrane from Run 3 is shown in Figure [H The 
backbone skeleton, plotted as lines, creates a hydrophobic network with chaotically scat- 
tered pores. These pores incorporate micelle-like clusters of sidechain sulfonates (shown as 
spheres), which are filled with water molecules and protons. The number densities p(r) of 
the hydrophobic part of the ionomer and of the absorbed water, averaged over a 100 fs run, 
are shown in Figure [2] and Figure [3l respectively. The density p(r) corresponds to the prob- 
ability of finding a particular membrane component, hydrophobic monomer or hydrophilic 
water, at the point r during a short simulation run. The quasi-regular network of polymer 
skeleton with interconnected hydrophilic pores changes its form slowly with time. 

The structure of the sulfonate multiplets was probed through the calculation of the 
sulfonate-sulfonate pair correlation function, 

V dns{r) , 
9 



Here dnsir) is the number of sulfurs located at the distance r in a shell of a thickness 
dr from a fixed sulfur atom. The function gss{r) indicates the probability of finding two 
sulfonate monomers at a separation distance r averaged over the equilibrium trajectory 
of the simulated system. Simulation results for (755 (r) for Runs 1-4 from Table [Til are 
shown in Figure [H The dry multiplets have no detectable internal structure except the 
strong maximum at r ^lAa. In hydrated membranes the correlation function gss shows 
shell-like oscillations, a recognizable fingerprint of solvation shells. The first maximum of 
gss{i^) corresponds to the closest-approach configuration between neighboring sulfonates. 
The second peak of gss{f^) stems from a configuration where two neighboring sulfonates are 
separated by single proton or water molecule. Finally, the third peak of gss{r) is related 
to configurations with more than one proton or water molecule between sulfonates. The 
dependence of the intensity of the correlations between the head groups on water content 
A is due to the dielectric screening properties of water: the more the water content in 
the membrane, the weaker the sulfonate-sulfonate interactions. Similar results have been 



reported in the simulation results of Refs. 
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51 



52| . The greatly reduced intensity of 



the first peak of (755 (r) at A = 5 can be understood as the onset of the improbability of the 
closest-approach sulfonate-sulfonate configurations in hydrated multiplets. 

We determine the size of a multiplet as the position of the global minimum (also 
known as the radius of the first coordination sphere) of the pair correlation functions gss 
in Figure [H This position depends on the membrane hydration level A, and can be used 
to calculate the number of head groups xs inside the multiplet according to the following 
relation 



Xs = i gss{r)^TTr dr. (9) 







The calculated values for the parameters and xs are given in Table UTTl for the membranes 
rji and 772- There is a clear indication of the fact that the multiplets shrink in size at the 
onset of membrane solvation, which corresponds to the transition from Run 1 to Run 2. This 
shrinking, which is not in accord with the classical theories of cluster swelling in ionomers, 
is accompanied by a multiplet splitting into smaller parts. For instance, the dry multiplet 
in the membrane 772 has a size R^ = 5a and consists of xs = 22 head groups. Following 
hydration by a water content as low as A = 1, this multiplet effectively splits into two 
smaller parts of size R^ = 4.2(T, each of them consisting of only 13 head groups. These 
small multiplets will consequently swell, keeping the number of their sulfonate population 
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constant, when additional water is absorbed into the membrane. The swelhng radius is 
largely determined by the competition between two different internal energies, the elastic 
energy of the backbone material and the electrostatic energy of the pendant groups. 

The pair correlation functions gss{^) for two different sulfonate concentrations rji and 
r/2 are plotted in Figure El When the parameter r] decreases, (as seen by a comparison of 
the thin and thick lines in Figure [5|), the intensity of sulfur-sulfur correlations increases. 
This effect stems from the interplay between the electrostatic screening length Id and the 
average separation distance I between the sulfonates. For an ionomer with a high sulfonate 
concentration rj, one generally has Id <l, and thus the electrostatic correlations between the 
sulfonates are negligible. In this case a nanophase separation in the membrane is possible 
only due to the hydrophobic/hydrophilic immiscibility between the backbone and sidechain 
segments of the membrane. In the opposite case, when rj is small and Id >l, the Coulomb 
correlations become sufficiently strong to force the sulfonates to form compact multiplets. 



B. Separation distance between multiplets 



It is a well established fact that the nearest-neighbor separation distance between the mul- 
tiplets and the connectivity of multiplets into a network of hydrophilic pathways are the main 
contributing factors to the transport properties of ionomers. The typical multiplet-multiplet 
nearest-neighbor distances can be directly deduced from the density-density correlations in 
the network of head groups by consideration of the structure factor. 



S{q) 
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J2 cos (qfl) 
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sin (qfl) 


1=1 




i=l 



(10) 



When there is no preferential ordering of the hydrophilic domains in the membrane, the 
structure factor of the sulfonates is isotropic, and hence depends only on the modulus q= \q\ 
of the wave-vector. The calculated structure factors S{q) for the dry and hydrated mem- 
branes are presented in Figure [6l The ionomer-peak position in the low g-region corresponds 
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40 



53 



Ml 



to the length of the density-density correlations R = 2n/q of sulfonates 

The nearest-neighbor distance between the multiplets can also be deduced, though less 
precisely, from the position of the long-range maximum i?£ of the pair correlation functions 
gss{r) in Figured! As in the multiplet splitting effect, corresponding to the reduction in 
the multiplet size at the onset of membrane solvation, the nearest-neighbor distance 



11 



also decreases to smaller values according to the results of Run 1 and Run 2. This is 
a consequence of the increase in the multiplet population ^ = Ns/xs- For example, the 
number of multiplets in the membrane rji increases from ^=30 in the dry membrane to ^=50 
in the hydrated membrane. The calculated values for R and -Rf, seen in Table UTTl match 
each other perfectly. We note that the increase of the average multiplet separation distance 
R for Runs 2-4 is a clear sign of membrane swelling, which is in accord with the results of 
Refs. Q,Q,Q- 



C. Swollen multiplets inside a water shell 

The pair correlation and the structure factor analysis, implemented in the previous sub- 
section, can be also exploited to examine the water clustering features in hydrated mem- 
branes for Runs 2-4. We calculate the size of the water cluster R^ from the water-water 
correlation gww{r) shown in Figure [71 The nearest-neighbor cluster separation distances 
R^ were evaluated from the water-water structure factors. The calculated values for both 
parameters are given in Table lllli There is good agreement between the water-water and 
the sulfonate-sulfonate multiplet nearest-neighbor distances R^ and i?f . This is an indirect 
verification of the fact that the ionomer cluster is a mixture of sulfonates and absorbed water 
molecules. The distribution of water molecules inside the ionomer cluster can be analyzed by 
comparing the water cluster size R^ with the sulfonate multiplet size Rj^. Whereas for the 
membrane rji there is an excellent match between these two parameters, for the membrane 
?72 the water clusters are consistently bigger than the sulfonate multiplets. Based on this 
result we conclude that a part of the total water loading per multiplet in fact exists outside 
the multiplet boundaries. This 'outer' water shell encapsulates the multiplet and facilitates 
the formation of narrow water channels between the swollen multiplets. These channels, 
clearly seen in Figure [Sj are the pathways through which the ionomer absorbs more solvent 
upon its hydration. The water channels are also an attractive place for the unclustered head 
groups, and assist the proton diffusion between neighboring multiplets. We remark that free 
bulk-like water would form in the interior of the multiplet only at sufficiently high solvation 
levels A |55i], a case not considered in this work. 

The average number of water molecules Xw per water cluster was calculated by using 
Eq.(l9l) for gwwif^)- This parameter, together with the parameter describing the water-per- 
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sulfonate ratio Xwlxs are given in Table UTTl For the membrane r]i we obtain Xw/xs = A, 
a predicted result for the ionomer cluster with i?^^ = R^. However for the membrane 772 
the ratio iw/is > A. This unexpected result can be interpreted in the following manner: 
when the sulfonate concentration approaches the percolation threshold for head-groups, a 
fraction of the sulfonates are randomly distributed between the existing multiplets. These 
bridging sulfonates cannot retain their full solvation shell with A water molecules in the 
hostile environment of hydrophobic backbones. The excess water molecules stripped from 
these 'bulk' sulfonates are consequently redistributed between the existing multiplets. This 
leads to the formation of an outer solvent shell around each multiplet. 

D. Swelling-induced ionomer deformation 

The polymer backbone and sidechains sustain conformational changes when the mem- 
brane swells. Two different types of deformation, an elongation (stretching) deformation 
and a coiling (frustration) deformation of polymer chains can be conveniently resolved using 
the probability distribution -P(a) of the dihedral angle along the polymer chains. 

The probability distribution P{a) of the dihedral angle along the sidechain is shown in 
Figure [8] for the membrane r/i. The sidechains have two gauche (±82 degrees) and one trans 
conformations. When the membrane absorbs water, the sidechains undergo a deformation 
in which a part of the gauche conformations transform into trans conformations. The overall 
effect of this structural deformation is a structural relaxation of the sidechains, perceived 
as a stretching - an impact schematically illustrated in Figure M We have also detected 
a similar stretching-like structural relaxation for the backbone segments. As seen from 
Figure [TOl the probabilities of the two gauche (±125 degrees) and single cis (±0 degrees) 
backbone conformations in the solvated membrane diminish when the hydration parameter 
A decreases. 

The extent of sidechain relaxation sensitively depends on the sulfonate concentration 77. 
In Figure [TT] we compare the P{a) curves for the two membranes. It is noticeable that the 
sidechains are in a more relaxed state in the membrane rji compared to the membrane r)2. 
This is a direct consequence of the fact that the smaller number of head groups inside the 
multiplet provide a more relaxed configuration for sidechains compared to the case when a 
larger number of sulfonates are immersed into a smaller multiplet. The dihedral frustration 
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of sidechains in the high-sulfonate membrane is schematically illustrated in Figure [121 A 
similar dihedral frustration has also been detected for the backbone polymer: in the high- 
sulfonate membrane the backbone segments adopt a more curly conformation. 



E. Proton Diffusion 

The proton mobility in solvated ionomers is strongly affected by proton-head group asso- 
ciation effects. On the one hand, this association localizes the protons near the head groups, 
and therefore decreases the rate of vehicular diffusion across the membrane. On the other 
hand, the localization effect increases the rate of the hopping diffusion of protons from sul- 
fonate to sulfonate. This so-called surface diffusion is believed to be additionally enhanced 
by the water-proton electrostatic interactions and the side-chain thermal fluctuations. The 
strength of the proton-sulfonate association is commonly evaluated in the terms of the pro- 
ton distribution around the head-groups. Our simulation results for the sulfur-proton pair 
correlation function gsuij') are plotted in Figure fT3l The first proton shell, seen as a very 
high peak on the left side of Figure[T3l originates from the attractive Coulomb forces between 
the protons and the SO3 groups. The second peak oi gsuir) on the right side of Figure [H 
arises from the proton shells of neighboring sulfonates in the multiplet. The condensation 
effect of protons on the sulfonates is noticeably stronger in the membrane rji than in the 
case of the membrane 772. As a consequence, the proton mobility in the membrane r]2 will 
be higher. 

The effect of a proton-sulfonate association also depends on the membrane hydration 
level A {5^: the association is weak for the hydrated membrane with A = 1, whereas it is 
strong for the membrane with A = 5. The position of the minimum of gsnif') corresponds 
to the position of the first maximum of gss{f) in Figured! 

The mean square displacements (msd) of protons for different membrane hydrations A 
are plotted in Figure [T4l As is expected from the proton delocalization effect in hydrated 



membranes, higher membrane hydrations result in larger proton displacements [49|. The 
msd result for the membrane rji is below the corresponding result for the membrane 772 for 
Run 4. This happens partly due to the strong proton delocalization effect, and partly due 
to the small nearest-neighbor distances in the membrane 772- 
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The calculated values for the diffusion coefficient of protons, 

msd(t) 



D 



lim 

t—>oo 



Qt 



(11) 



are gathered in Table UTTl The proton diffusion, similar to the proton mobility, is stronger in 
membrane 772 than in membrane rji because of the low proton-sulfonate association. There 
are two other factors that contribute to the proton diffusion of membrane 772: the existence 
of 'bulk' sulfonates between neighboring multiplets and the accumulation of water molecules 
around the multiplets. Both these factors can lead to the formation of temporary bridges, 
sulfonic and/or solvent in nature, between the multiplets. Our results for proton diffusion 



49|. However they are small com- 



are in good agreement with the simulation results of Ref. 
pared to the proton diffusion coefficients experimentally observed in fully hydrated Nafion 
ionomers. This discrepancy is most probably not due to our neglect of the Grotthuss mech- 



anism, which is strongly suppressed when A is small 
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581, but is a consequence of the 



reduced number of pathways in our low-humidity, low-sulfonate model. 



V. DISCUSSION 



We have investigated the swelling properties of multiplets in low humidity ionomers with 
low sulfonate concentration by examining different models for the sidechain architecture. 
Our primary goal was to determine the dependence of multiplet swelling on the hydration 
level A and the sulfonate concentration t] of the membrane. 

Our main result is the fission of the sulfonate multiplets into smaller parts at the onset 
of membrane hydration. This behavior is not explained by the classical theories of cluster 
swelling in ionomers, according to which the swelling should be a continuous and monotonic 
process of multiplet expansion. The resultant small multiplets will consequently swell, keep- 
ing the number of their sulfonate population constant, when more water is absorbed into the 
membrane. We have also found that the location of the of water in low-sulfonate membranes 
strongly depends on the sulfonate concentration. For a relatively low sulfonate concentra- 
tion nearly all sulfonate groups are in multiplet formations. The average water loading 
parameter per multiplet Xw/Xs, where xw is the number of water molecules belonging to 
the multiplet, and xs is the number of sulfonates in the multiplet, perfectly matches the 
water content of the membrane A. However, for relatively high sulfonate concentrations, the 
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water loading parameter per multiplet Xwlxs is consistently larger than the parameter A 
for the membrane hydration levels considered. We assume that, when the sulfonate concen- 
tration approaches the percolation threshold for head-groups, a fraction of the sulfonates 
are randomly distributed between the existing multiplets. These bridging sulfonates cannot 
retain their full solvation shell in the hostile environment of hydrophobic backbones. The 
excess water molecules stripped from these 'bulk' sulfonates are consequently redistributed 
between the existing multiplets. The results of our structural analysis confirm the formation 
of unexpected water shells around sulfonate multiplets. The multiplet fission and the water 
encapsulation effects are illustrated schematically in Figure fTSl 

Our discovery of the uneven distribution of the water-to-sulfonate loading in the ionomer 
opens a new window into the percolation characteristics of the hydrophilic network in 
ionomers. It is no longer sufficient to have a continuous pathway among sulfonates in 
order for percolation of protons to occur, as some of these sulfonates may be found in the 
hydrophobic material, where they are not capable of contributing to proton transport. The 
predicted hydration levels necessary for good transport of protons will thus be higher than 
they would be if the presence of sulfonates encapsulated in backbone material were ignored. 

We have also analyzed the structural deformations occurring in the ionomer as a result of 
membrane swelling, and found that in swollen membranes the ionomer is in a more relaxed 
state. The degree of relaxation, however, is sensitive to the sulfonate concentration: the 
sidechains and backbones are found to be more relaxed in low-sulfonate membranes. This 
result is a direct consequence of the fact that in low-sulfonate membranes the sulfonate 
cluster is less dense, and can relax more readily than in the denser environment of the 
high-sulfonate material. However, proton diffusion is stronger in high-sulfonate membranes, 
and can potentially benefit from the formation of temporary solvent and sulfonate bridges 
between the multiplets. 

An interesting question, yet to be resolved, is whether the multiplet splitting and shrink- 
ing effects depend on the pretreatment history of the membrane. The membrane morphology 
is known to be affected by the type of pretreatment, such as boiling, annealing, drying, pol- 
ing, stretching, etc., and by the order in which these steps are taken. In most cases the 
impact of the pretreatment is either the formation of a new morphology with an anisotropy 
in the backbone and sidechain orientations, or the reshaping of the network of hydrophilic 
clusters. In our current work the dry membrane was 'numerically pretreated' by our frag- 
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mentation and de-fragmentation procedures, as described in section III. We assume that our 
membrane has a network of hydrophihc pores resembhng the network in a mold-extruded 
membrane, provided it has then been annealed. 

In order to analyze the consequences of the residual anisotropy in the ionomer on the 
multiplet reorganization effects reported in this work, we carried out test simulations for a 



poled and dry Nafion-like ionomer. According to our previous results on ionomer poling [33|, 
rod-like aggregations of head groups are formed along the direction of the applied electric 
field. The poled structures were found to be stable after the release of the poling field. 



One of the poled structures of Ref. [33| was used as a starting configuration for Run 1 of 
our current work. Our simulation result indicated that a similar reorganization effect of 
sulfonate multiplets, as seen in the case of isotropic membranes, takes place. Hence, we 
conclude that the splitting and shrinking effects are robust against structural anisotropy in 
the membrane. 

We also performed test simulations to clarify the nature of multiplet reorganization in 
dry membranes that had been previously swollen. The hydrated membrane from Run 4 
with water content A=5 was first dried through a simple elimination of all water molecules 
in the simulation box. Then the water-free membrane was gradually shrunk to the system 
size used for Run 1. The results obtained show that the initially dry membrane, membrane 
I, and the pretreated dry membrane, membrane II, have different structures. In the latter 
the multiplet sizes were smaller and close to the multiplet sizes corresponding to Run 2. 
However, after annealing at high temperatures, the discrepancies between the membranes 
disappeared, and both membranes exhibited the splitting and shrinking of multiplets at the 
onset of hydration. 

In future work, we plan to extend the model presented here to take into account the 
partial charges on the side-chain monomers. Our preliminary results indicate that a partial 
delocalization of the negative charge along the sidechain head group has a noticeable role in 
the membrane swelling process. 
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Table I: List of key variables. 

A number of water molecules per sulfonate group 

a monomer diameter 

e dielectric constant of medium 

Qp,Qs,Qo3 normalized charges of protons, sulfur atoms and head group oxygens 

El J Lennard- Jones interaction parameter between monomers 

fcfe, kd, ka stretching, bending and torsion force constants 

ro equilibrium bond length for backbone and sidechain 

^0 equilibrium bending angle for backbone and sidechain 

kB, T Boltzmann constant and system temperature 

n\,n2 number of hydrophobic and hydrophilic monomers per sidechain 

713 number of hydrophobic backbone monomers between adjacent sidechains 

ry molar concentration of head groups 

Nq Avogadro's number 

N , Ng total number of ionomer monomers, number of sulfonates 

L, V length of simulation box, volume of simulation box 

gssi^), gsnir) sulfonate- sulfonate and sulfonate-proton pair correlation functions 

gwwi'f) water- water pair correlation function 

-^A! ^i^^ of sulfonate-sulfonate multiplets and water clusters 

R%, nearest-neighbor distance between for sulfonate multiplets and water clusters 

XSi Xw number of head groups and number of water molecules in a multiplet 

^ number of multiplets in the simulation box 

Id electrostatic screening length 

I average separation distance between sulfonates 

S{q), q = 2-K/r sulfonate-sulfonate structure factor 

R correlation length of density-density fluctuations of sulfonates 

Pia) Probability distribution of the dihedral angle along the polymer segments 

D diffusion coefficient of protons 

p 3D density of membrane components 
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Table II: Parameters used in simulation runs. Here ni + n2 is the total number of monomers per 
sidechain, A is the water content per head group. 



Runs hydration model ni + n2 

Run 1 dry ionomer with no water, A=0 10 

Run 2 ionomer with explicit water, A=l 9 

Run 3 ionomer with explicit water, A=3 9 

Run 4 ionomer with explicit water, A=5 9 
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Table III: Calculated ionomer parameters for membranes with sulfonic molar concentration rji and 
r/2 (shortly called membranes rji and r]2 in the text. The definitions of the parameters used are 
given in Table [B 



membrane rji 



A 




i?| 


R 




Rq 


XS 


XW 


XwlXS 


Z)(cm^/sec) 





5.3 


9.6 


10 
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1 


4.7 


7.9 


7.5 


4.7 


7.9 


10 


10 


1 


2.7x10-6 


3 


4.9 


8.1 


7.7 


4.9 


8.1 


10 


30 


3 


4.5x10-6 


5 


5.0 


8.5 


8.0 


5.0 


8.5 


10 


50 


5 


5.2x10-6 



membrane 772 



A 


Ri 




R 




R^ 


XS 


XW 


XwlXS 


Z)(cm^/sec) 





5.0 


9.4 


9.97 
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1 


4.2 


7.9 


7.64 


4.8 


7.7 


13 


17 


1.2 


3.5x10-6 


3 


4.4 


8.0 


8.0 


5.1 


8.04 


13 


62 


4.5 


5.2x10-6 


5 


4.6 


8.2 


8.35 


5.2 


8.32 


13 


101 


7.8 


6.5x10-6 
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110 

Figure 1: (Color online) A typical snapshot of hydrated membrane 772 from Run 3. The spheres 
represent the end-group oxygens of the sidechains. The polymer is shown by red lines. Different 
bead colors correspond to different bead altitudes, with a blue color for low-altitude beads (at the 
bottom of simulation box) and a red color for high-altitude beads (at the top of simulation box). 
The size of all structural elements is schematic rather than space filling. The water molecules and 
protons are not shown for the sake of clarity. 
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Figure 2: (Color online) 3D density p{'r) of the hydrophobic part of the membrane r]2 for Run 3. 
The color gradient from dark blue (black in printed version) to dark red (gray in printed version), 
corresponds to the variation of membrane density from low to high values. The axis dimension is 
in A. 
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Figure 3: (Color online) 3D density of water channels for the membrane 772 and Run 3. The color 
gradient from dark blue (black in printed version) to dark red (gray in printed version) corresponds 
to the variation of water density from low values to high values. The axis dimension is in A. 
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Figure 4: (Color online) Sulfonate-sulfonate pair correlation function 55s(r) for membrane 772 as a 
function of sulfur-sulfur separation distance r for Runs 1-4 from Table [III Solid line with circles 
- Run 1, solid lines with squares- Run 2, dashed line- Run 3, full line- Run 4. The bottom figure 
shows in detail the long-range tail of gss{f ) used to determine the average multiplet size and the 
separation distance between the sulfonate multiplets R^. The calculated values for the parameters 
R^ and R^ are given in Table IIIII 
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Figure 5: (Color online) Sulfonate-sulfonate pair correlation function gss{i") for the membranes 771 
and r/2 as a function of the sulfur-sulfur separation distance r. Lines with symbols- Run 1, full 
lines- Run 4. 



29 




Figure 6: (Color online) Small angle ionomer peak region of sulfonate-sulfonate structure factor 
S{q) for the membranes r/i and r]2. Line with circles- Run 1, line with squares- Run 3, full line- 
Run 4. The step size 6q = 2-k/L « O.OSA"^ defines the resolution along the z-axis. 
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Figure 7: (Color online) Water- water pair correlation function gww{i^) for the membrane 772 as a 
function of the water- water separation distance r for Runs 2-4 from Table HTl Solid line with circles 
- Run 2, solid lines with squares- Run 3, full line- Run 4. The arrows show the average multiplets 
size B}^ and the separation distance between the water clusters. The calculated values for the 
parameters and are given in Table lllli 
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a (degrees) 



Figure 8: (Color online) The probability distribution P{a) of the dihedral angle along the sidechain 
of membrane rji, and for Run 2 and Run 4. The areas of the gauche conformation (a = 82 degrees) 
and the trans conformation (a = 180 degrees) are shown separately. 
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Figure 9: (Color online) Schematic pictures explaining the sidechain stretching-like relaxation as a 
result of multiplet swelling from A = 1 to A = 5. The small hollow spheres are the water molecules, 
gray (yellow in online version) small spheres with attached tails are for sidechains, big spheres 
represent the multiplets. 
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a [degrees] 

Figure 10: (Color online) The probability distribution P{a) of the dihedral angle along the backbone 
of membrane rji for Run 1 and Run 4. The areas of a cis conformation (a = degrees) and the 
gauche conformation (a = 125 degrees) are shown separately. 
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Figure 11: (Color online) The probability distribution P{a) of the dihedral angle along the 
sidechains of membranes r/i and rj2, and for Run 4. The areas of the gauche conformation (a = 82 
degrees) and the trans conformation (a = 180 degrees) are shown separately. 
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Figure 12: (Color online) Schematic pictures explaining the differences between the stretching-like 
relaxation of sidechains for the membranes rji and r]2. The small hollow spheres are the water 
molecules, gray (yellow in online version) small spheres with attached tails are for sidechains, big 
spheres represent the multiplets. 



36 




Figure 13: (Color online) Sulfonate-proton correlation function gsH{f) as a function of the sepa- 
ration distance r. The first and second peak areas are shown separately. Thick lines are for the 
membrane r^i, thin lines are for the membrane r/2. Line with symbols- Run 3, full lines- Run 4. 
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Figure 14: (Color online) Mean squared displacement of protons as a function of time for Runs 2-4. 
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Sulfonate concentration 



Figure 15: Schematic illustration of the multiplet hydration. At a low sulfonate concentration the 
hydrated multiplets consist of sulfonates and water molecules. At a higher sulfonate concentration 
each of the multiplets is surrounded by a water shell. The splitting of dry multiplets into smaller 
hydrated multiplets is also sketched. Vertical/horizontal hatching is used for the water (W) and 
the sulfonate (S) areas of the multiplet. 
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